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Abstract 

We consider the xo-dependence of the density of states inside the normal 
metal of a superconductor - normal metal - superconductor (SNS) junction. 
Here xo is the phase difference of two superconductors of the junction. It is 
shown that in the absence of electron-electron interaction the energy depen- 
dence of the density of states has a gap which decreases as xo increases and 
closes at xo = ^- Both the analytical expressions for the xo dependence of 
the density of states and the results of numerical simulations are presented. 
PACS index category: 05.20-y, 82.20-w 

I. INTRODUCTION 

The proximity effect in a normal metal near superconductor-normal metal boundary has 
been studied for many years. Recently, progress in microfabrication technology has revived 
the interest in this phenomenon. Due to Andreev reflections at superconductor-normal metal 
boundaries, the density of states inside the normal metal of a superconductor-normal metal 
(SN) junction is different from that of the bulk normal metal. This phenomenon has been 
studied both theoretically and experiment ally^Hi. In the superconductor - normal metal - 
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superconductor (SNS) junction it has been shown that at xo = the Andreev reflections 
at the superconductor-normal metal interfaces cause the density of states to have a gap, 
which is of the order of the Thouless energy E c = D/L 2 <C Ao- Here xo is the phase 
difference between the two superconductors, Ao is the modulus of the order parameter in 
the superconductors, L is the length of the normal metal of the junction (see Fig.l) and D 
is the diffusion constant in the metal. 

In this paper we study the xo dependence of the density of states z/(e, xo) m the case 
when the length of the normal metal region is much larger than the electron elastic mean 
free path (diffusive regime). We study both analytically and numerically the behaviour of 
the density of states near the gap edge. This quantity can be measured for example in a 
tunneling experiment. We will show below that it also determines the conductance of the 
junction Gsns a t low temperatures T. 

The density of states in the normal metal can be expressed in term of the retarded 
Green's function g R (e,x) = cos#(e, x)i 



u r L / 2 

v(e) — — I cos#i(e, x) cosh^e, x)dx. (1) 



Here vq = mpp is the density of states in the bulk normal metal and pp is the Fermi 
wave length; 6 = 6i + i9 2 and x — Xi + % are the complex variables described by Usadel 
equationsi, 

—d x 2 6(e, x) + ie sin 9(e, x) — —(d x x{e, x)) 2 sin 2#(e, x) = A N (x) cos 9, 

d x {d xX {t,x)sm 2 6{e,x)} = 0, (2) 

Aat(x) = 7tv J decos0i(e, x) sinh^e, x) tanh(e/2/cT). (3) 

Aat and 7tv > are the modulus of the order parameter and the dimensionless repulsive 
interaction constant inside the normal metal. 

The boundary conditions for Eq.2 are (we consider the case when the transmission coef- 
ficient of the normal metal-superconductor boundary t — 1), 



X(e,x = ±-) = ±y. 



(4) 



Consider the situation when there is no electron-electron interaction in the normal metal 
and A N = 0. We will see that in this case the energy gap E g (xo) of the electron spectrum is 
a decreasing function of xo and is equal to zero only at xo = ^ (See Fig. 2). In two limiting 
cases, when xo is close to or it we have 



E g ( Xo ) = E c 



C 2 (l - C lX l) Xo < vr 
Csiir - xo) 7T - xo < 7T 



(5) 



II. THE ENERGY GAP 



To get this result we take into account that at e < E g (xo), #i = tt/2 following Eq.l. In 
this case the second integral of Eq.2 gives the solution for 2 (e, x): 

r02o f i 2t 

/ d# 2 { — (sinh 6 20 - sinh # 2 ) + a 2 (e) (cosh" 2 6 2 - cosh" 2 # 20 ) }" 2 = ( 6 ) 

•/0 2 (e,x) E c L 

where 6* 20 (e) and a(e) are the functions of energy e and Xo determined by the equations 

f 920 d6 2 {-|- (sinh 6 20 - sinh 9 2 ) + ^^(cosh' 2 # 2 - cosh" 2 # 20 )}~^ = 1, 

a{e) r 20 ^ 2 cosh- 2 ^ 2 {^(sinh^ 2 o-sinh^ 2 ) + ^^(cosh- 2 ^ 2 -cosh- 2 ^ 20 )}-^ = ^. 
Jo E c 2 2 

(7) 

One can show that Eq.7 has solutions only at low energy e. The energy gap E g (xo) 
corresponds to the maximum value of e, at which a solution exists. Beyond this value, 
0\ 7^ 7r/2 and therefore the density of states z/(e) is non zero. 

Let us first consider the limiting case corresponding to xo <^ tt- We can expand Eq.7 
with respect to xo an d as a result we have the following equation: 
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e 



7 = r° d9 2 (smh6 20 - sinh0 2 )-* - (^-fA(9 20 ), (8) 
c Jo 2 



where A is of the order of unity. The right hand side of Eq.8 as a function of 9 2 o has 
a maxima equal to y/U^il — Cixjj/2). Here C\ >2 are numerical factors of order of unity. 
Their values can be obtained from the numerical solution of the Eq.20: C\ = 0.91 and 
C 2 = 3.122. At low energies when e <C E c , Eq.8 has a solution 9 2 q ~ t/E c while at 
large energies e > C 2 E C (1 — CiXo)> Eq.8 does not have a solution. Therefore in this limit 
E g = C 2 E C {1 — CiXq), as given in Eq.5. 

When xo — 71 ^ tt> the energy gap is small, i? s <C _E C and we can expand Eq.7 with 
respect to the small parameter e/E c . As a result we have 

7r 2 tanh02o{8£i + (tt - Xo)^ sinh0 2 o} _1 - (9) 



(vr - xo)-#c 

where -B1.2 are also constants of the order of unity. Again the right hand side of Eq.9 as a 
function of 9 2 q reaches a maxima when 

cosh#20 = B (n - Xo)~^ (10) 

Therefore, at e > CsE^tt — xo), Eq.9 has no solution and E g (xo) is given by Eq.5. From 
the numerical solution we find Bq = 2.21 and C3 = 2.43. The insert of Fig. 2 shows the 
linear dependence of the energy gap near Xo = n - 

III. THE GAP EDGE 



Let us turn to the calculations of the e-dependence of the density of states u(e, xo) at 
e>E 9 . 

In the region e — E g <^ E g , the quantities 

59 1 (e,x) = 0i(e > E g ,x) - |, 
59 2 (e,x) = 9 2 (e> E g ,x)-9 2 (E g ,x), 
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<*Xi(e,aO = Xi(e > - Xi(E g ,x), 

Sx2(e,x) = X2(e > E g ,x), (11) 

are small and go to zero as e approaches E g . Here 9 2 (E g ,x), Xi(E g ,x) are the solution of 
Eq.2 at e = E g given by Eqs.6,7. When n — xo <C 7r 

sinh 9 2 (E g , a;) = sinh # 20 cos — , 

^Xi(g g ,x)= ""f^ r , (12) 
A v 9 ' ; 2Lcosh 2 fl 2 (£ 9 ,a;) V ; 

where 9 2 o is given in Eq.10. 

Expanding Eq.2 with respect to 59 ij2 , 8x1,2 we the following set of equations 

{jd 2 x + E g smh9 2 + j(d xX i) 2 cosh29 2 }56 1 - ^d X i sinh 2^X2 

-S%(^smh9 2 + ±D(d xX i) 2 cosh26 2 ) 

+59 1 {E g cosh 9 2 S9 2 + D(d xX i) 2 sinh 29 2 59 2 + DdJxAxi cosh2fl 2 - j(d x 5 X2 ) 2 cosh 29 2 ) 

+5 2 6»iD9 :E (5x2<9 a; Xi sinh 26*2 - DdJx259 2 d x x\ cosh26> 2 - ^d^Xi^W sinh20 2 

= -(e-e p )sinh0 2 <S0i (13) 



{^d 2 x + E g sinh9 2 + j(d xX i) 2 cosh29 2 }59 2 + ^d xX i smh26 2 d x 6 X i 
-^(^coshfc + j(d xX i) 2 sinh 2fl 2 ) + 89 1 DdJ X 2d x Xi cosh20 2 

-j(^5% 2 ) 2 sinh 2^ 2 = -(e-e 9 )coshfl 2 (14) 



0{dJx2 cosh 2 # 2 - £0id*Xi sinh 2^ 2 } = (15) 



d x {dJxi cosh 2 9 2 + 59 2 d x xi sinh20 2 - 59\d x xi cosh20 2 + 56x8x5x2 sinh20 2 } = (16) 
with the boundary conditions for 59± t2 ,5xi,2' 

86 1 (e,x = ±|) = 59 2 (e,x = ±^) = 0, 

5 X i(e,x = ±^) = 5 X 2(e,x = ±^) = 0. (17) 



Consider the linear part of Eq. 13-17 at e — E g = 0. Its solution determines the spatial 
dependence of 59^2, 5xi,2- Averaging Eq.13-17 over the sample and comparing the nonlinear 
terms in Eq.13 with right hand side term proportional to (e—E g ) sinh#2<5#i, we can determine 
the value of 59 \ as a function of e — E g . As a result we have 



where C4 is a function of E g and consequently of xo- Substituting Eqs.10, 18 into Eq.l we 
obtain the asymptotic form valid when e — E g <C E g : 



with C4 ~ (n — Xo) _/3 ; where (3 ~ 2/3. The value of (3 is a result of a numerical calculation 
over the whole energy rangeB. The square root energy dependence of the density of states is 
illustrated in Fig. 3 which shows the square of the density of states as function of the energy 
very close to the energy gap. At higher energy (Fig. 3 insert), we find that the density of 
states exhibits a smooth bump above the gap when phase difference is small. This smooth 
maximum turns to a sharper and sharper peak as the phase difference approaches n, i.e. as 
the gap closes. 

It should be noted at this point that the asymptotic density of states at high energy 
or at Xo = ^ is smaller than the normal state value. The region near the superconductor 
provides a small contribution to the (spatially averaged) density of states. This apparent 
deficit of states is balanced by the excess density of states above the superconducting gap 
A which according to our assumption is far above the energy range of interest. This deficit 
does not exist at the center of the normal metal or in a geometry where the relative area of 
the boundary between N and S goes to zero. The latter case is met in the billard geometry 




(18) 



v{e, Xo) = ^oC^Xo) 



e-E g ( X o) 



(19) 



of Rei. 



IV. DISCUSSION 
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Although the above results were derived for the junction geometry shown in Fig.l, the 
Xo an d e dependences of the density of states in Eqs.5,19 as well as the statement that the 
gap closes at Xo = * are general. We believe that they are independent of the transmission 
coefficient of the superconductor-normal metal boundaries and the geometry of the normal 
region. The values of E c and Cj, however, depend on these parameters. Similar conclusions 
were also reached in@. 

So far we assumed the electrons inside the normal metal do not interact with each other. 
In the presence of electron-electron interactions, An is not zero and one can study the 
interaction effects on the density of state by taking into account An in Eq.7 for E g . Since 
7at <C 1, one can carry out the perturbative calculation with respect to jn- As a result, at 
Xo 7T, the gap turns out to be smaller than that given in Eq.5 for the noninteracting case, 
i.e. E g (0) — E 9 (jn) ~ 1nE c . Furthermore, the gap closes at x* smaller than n 



Quantum fluctuations of the phase of the order parameter can also change the results derived 
above. 

At finite temperature, in principle one also has to take into account the electron level 
broadening due to inelastic scattering. Such an inelastic process will introduce a temperature 
dependent density of state at the Fermi surface, i.e. v(e = 0) ~ 1/E c r t . However at T <C 
E g (xo), the inelastic scattering rate r" 1 due to electron-phonon interaction is exponentially 
small. 

Finally, let us calculate the conductance of the junction at T <C E c . The conductance 
Gsns is the proportionality coefficient between the applied voltage and the disspative current 
averaged over the period of Josephson oscillations. At low temperature, the main contri- 
bution to Gsns comes from the Debye relaxation mechanismll. The qualitative picture is 
the following. When voltage V is applied across the junction, the time dependence of xo is 
determined by the Josephson relation 



7T - X* ~ In 



(20) 



dxo 



2eV 



(21) 



dt 
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At small V the time dependence of v is determined by the corresponding time dependence 
of Xo{t)- I n other words the quasiparticle energy levels move adiabatically with frequency 
2eV. The electron population of the energy levels follows the motion of the levels and as a 
result the electron distribution becomes nonequilibrium. The relaxation of the nonequilib- 
rium distribution due to inelastic processes leads to the entropy production and therefore 
contributes to the conductance. As a result, we have!, 

Gsns = «tanh-|- d Xo r e ( Xo ){ f de '^2*Ay (22 ) 

Ti vq J Zkl Jo J-00 axo 

where r e (xo) is the energy relaxation time of quasi particle of energy e. In this case0, 

— ^— ~ jL exp (- Eg ( X o)/kT) (23) 

T e{X0) * l D 

is exponentially small in the time interval when E g (xo{t)) 3> T. is the Debye frequency. 
On the other hand, the concentration of quasiparticles in this case is also exponentially 
small. These two exponential factors cancel each other and the main contribution to Eq.22 
comes from the time interval when E g ( X o{t)) ~ T. Since Xo changes linearly with time and 
Eg(xo) vanishes linearly as a function of xo — tt at xo close to it, following Eq.22 we have 

Gsns « GVC 4 2 (^)r m T (24) 

Here G N = e 2 Du f and t^ 1 ~ T 3 /tt 2 D . C A is given in Eq.19. At E c ~ T, Eq.24 matches 
the result of the conductance, G^Ti n E 2 /T ', which was obtained mi at the high temperature 
limit T > E c . 

Following the above arguments, the time dependence of the conductance has the form of 
narrow peaks with the amplitude of order of GatC|(T / E c )r in E c and duration of the order of 
(eV)~ 1 TE c . This phenomenon can be connected with the well known cosxo problem which 
has been investigated both experimentally and theoreticallj^ 1 1 



V. CONCLUSION 
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In conclusion we have shown that the density of states in the normal part of a SNS 
junction has a gap which closes when the superconducting phase difference is Xo = 7r - 
The energy dependence of the density of states near the gap edge has been calculated. We 
should mention that nonequilibrium effects in superconductors which are connected with 
time dependence of quasiparticle spectrum have been considered in0'@. In these papers the 
enhancement of the critical current due to the nonequilibrium effects was considered. The 
effect considered there is proportional to (eW e ) 2 while the DC current calculated above is 
linear in (eVr e ). 

We would like to acknowledge discussions with H. Courtois. This work was partially 
supported by the NATO CRG 960597. 
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FIG. 1. The S-N-S system: the origin of coordinate x is the center of the normal metal (length 
L). The phases of the superconducting order parameters in the right and left superconducting 
electrodes are respectively — Xo/2 and +xo/2. 
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FIG. 2. Energy gap E g in units of the Thouless energy E c vs the phase difference xo between 
superconducting contacts. The insert shows the linear dependence of E g near xo = 
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FIG. 3. The square power of the density of states vs reduced energy e/E g (xo) near the gap 
edge. Curves labeled 0, 1, 2, 3, 4, 5, 6 and 7 are for xo = 0, tt/2, 37r/4, 7it/8, 157r/16, 3l7r/32, 
637r/64 and it. The insert shows the full density of states curves (here E g <C Aq). 
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FIG. 4. Divergence of the C4 numerical coefficient near xo = The dashed line is an 
asymptotic law (it — xo)~^ with (3 = 2/3. 
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